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Abstract 

Using time-resolved photoluminescence spectroscopy over a wide range of temperatures, we were able to 
probe both radiative and nonradiative relaxation processes in luminescent porous silicon. By comparing the 
photoluminescence decay times from freshly prepared and oxidized porous silicon, we show that radiative processes 
should be linked with quantum confinement in small Si nanocrystallites and are not affected by oxidation. In contrast, 
nonradiative relaxation processes are associated with the state of oxidation where slower relaxation times characterize 
hydrogen-terminated porous silicon. These results are in a good agreement with the extended vibron model for small 
Si nanocrystallites. 
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Background 

The efficient room-temperature visible photolumines- 
cence (PL) from porous silicon (PSi) has attracted much 
attention in recent years, mainly due to open questions 
and controversies concerning the mechanism responsible 
for the PL emission [1-7]. In addition, numerous PSi-based 
devices having potential applications in diverse fields such 
as photonics, optoelectronics, and photovoltaics, were 
proposed and investigated [8-15]. In particular, PSi has 
been considered as an attractive candidate for sensing ap- 
plications [16-21] where its large surface area can be 
exploited for enhancing the sensitivity to surface interac- 
tions. In such a sensor, the PL emitted from PSi can be 
used as a transducer that converts the chemical inter- 
action into a measurable optical signal. For example, PL 
quenching due to surface interactions with various chem- 
ical species has been utilized for developing various bio- 
photonic sensors [16,22,23]. 

Originally, the efficient PL from PSi was attributed 
to quantum confinement (QC) of charged carriers in Si 
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nanocrystallites located in the PSi matrix [24]. Experi- 
mental evidences supporting this model include a shift of 
the energy bandgap with size [1-3,25,26], resonant PL at 
low temperatures [27-29], and PL decay time spectroscopy 
[1,2,27]. However, the QC model cannot account for other 
experimental observations, mainly the dependence of the 
PL on surface treatments [30-34]. Several reports pro- 
posed a more complex picture of QC combined with 
localization of charged carriers at the surface of the nano- 
crystals [35-38], particularly the work of Wolkin et al. [36] 
who demonstrated a strong dependence of the PL on sur- 
face chemistry. This group has shown that while in fresh 
PSi the PL peak energy depends on the size of the nano- 
crystals (i.e., follows the QC model), the QC model cannot 
account for the limited PL shift observed for oxidized PSi. 
By introducing surface traps into the model, the behavior 
of the PL peak energy for oxidized PSi could be explained 
[36]. Other reports have shown that both QC and surface 
chemistry shape the PL characteristics [37,38]. The ex- 
tended vibron (EV) model provides a simple explanation 
to the mutual role of surface chemistry and QC [39-41], 
According to this model, QC affects radiative processes 
that are less sensitive to the state of the surface, while 
nonradiative relaxation processes are mostly influenced by 
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the surface chemistry. However, both QC and surface 
chemistry contribute to the efficient PL from PSi. 

In this work, we investigate the role of surface chemis- 
try, particularly the relationship between the state of oxi- 
dation and the PL characteristics of luminescent PSi 
samples. We examine the contribution of radiative and 
nonradiative decay processes to the overall PL lifetime 
and the sensitivity of these processes to surface treat- 
ments. Furthermore, we examine the EV model by com- 
paring radiative and nonradiative decay times of freshly 
prepared hydrogen-terminated PSi (H-PSi), with those 
of oxidized PSi (O-PSi). This allows to experimentally 
test the hypothesis that radiative processes are not sensi- 
tive to surface treatments while nonradiative processes 
are. Utilizing temperature-dependent, time-resolved PL 
(TR-PL) spectroscopy [42], we extend our previous work 
on silicon nanocrystals embedded in Si0 2 matrices and 
silicon nanowires [37,41,43,44] to PSi, as this system 
allows a modification of the surface chemistry by sim- 
ple means and tracing quite accurately the state of the 
surface. 

Methods 

PSi samples were prepared by electrochemical etching of 
p-type (10 to 30 H em) silicon wafers under standard 
dark anodization conditions [25,26]. A 1:1 mixed solu- 
tion of aqueous hydrofluoric (HF) acid (49%) and etha- 
nol was used as the electrolyte at a current density of 
70 mA cm" 2 for 200 s to yield a PSi layer of approxi- 
mately 9.5 um (measured by scanning electron micro- 
scope) with average pore size of a few nanometers [25]. 
The freshly prepared PSi is terminated by Si-hydrogen 
bonds that are known to be quite unstable under ambi- 
ent conditions. These bonds are subsequently replaced 
by the more stable Si-oxygen bonds upon exposure to 
air. Hence, in order to investigate the optical properties 
of H-PSi, we introduced the freshly prepared samples 
into a vacuum optical cryostat and kept them under 
vacuum conditions for the entire experiment. Oxygen- 
terminated O-PSi was obtained after taking the same PSi 
sample out of the vacuum cryostat and letting it age under 
ambient conditions for 6 days. The state of the PSi surface 
(having either Si-O or Si-H bonds) was monitored by Fou- 
rier transform infrared (FTIR) spectroscopy. To eliminate 
interference phenomena, thinner PSi samples were pre- 
pared for these measurements (10 s of anodization under 
the same conditions, resulting in approximately 450 nm 
thick PSi film). Bruker's Vertex- V70 vacuum FTIR spec- 
trometer (Bruker Optik GmbH, Ettlingen, Germany), 
equipped with a mercury-cadmium-telluride (MCT) pho- 
tovoltaic detector, has been exploited for these experi- 
ments. Measurements were performed in the grazing 
angle reflection mode, at an incidence angle of 65° and 
under p-polarization (to enhance the sensitivity to surface 



bonds [45]). For continuous wave (cw) PL and TR-PL 
measurements, the samples were excited by Ar + ion laser 
operating at 488 nm while the PL signal was dispersed by 
a 1/4-m monochromator and detected by a photomulti- 
plier tube. For time-resolved measurements, the laser 
beam was modulated by an acousto-optical modulator 
driven by a fast pulse generator, while the PL signal has 
been analyzed by a gated photon counting system. During 
PL measurements, the samples were kept under vacuum, 
in a continuous-flow liquid helium optical cryostat that al- 
lows temperature control from approximately 6 K up to 
room temperature. 

Results 

IR absorbance spectra of H-PSi (red line) and O-PSi 
(black line) are presented in Figure 1. Si-OH and Si-O-Si 
vibrational bands at 875 cm" 1 and 1,065 to 1,150 cm -1 
respectively [46-48], which indicate the presence of oxy- 
gen in the films, clearly increase after 6 days of exposure 
to ambient atmosphere. The Si-hydrogen vibrations 
[46,47] at 906 cm" 1 and 2,112 cm" 1 did not fade away 
after 6 days of exposure, and typically disappear only 
after several weeks of exposure to air [49]. Hence, these 
IR absorbance spectra confirm the modification of the 
PSi's surface during the exposure to air. 

The cw-PL spectrum of H-PSi, measured at room temp- 
erature with a PL maximum at approximately 1.80 eV 
(about 690 nm) and a full width at half maximum (FWHM) 
of about 0.4 eV, is presented at the inset to Figure 2. A 
similar spectrum with a slight blue shift of the PL max- 
imum to 1.85 eV (approximately 670 nm) has been mea- 
sured for O-PSi, in agreement with results obtained in 
references [50-52]. In order to probe both radiative and 
nonradiative relaxation processes, the PL decay curves 




Figure 1 FTIR spectra. Infrared absorption spectra of H-PSi (freshly 
prepared PSi) and O-PSi (the same sample after aging). Main Si-H, 
Si-OH, and Si-0 vibration modes are marked. 
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Figure 2 PL decay curves. The PL decay curves of H-PSi measured 
at a photon energy of 2.03 eV (610 nm) and at various temperatures. 
The solid lines present the best fit to a stretched exponential 
function (Equation 1). Inset shows the PL spectrum of H-PSi 
measured at room temperature. 



were measured at several photon energies and at tem- 
peratures ranging from 6 K up to room temperature. As 
will be discussed and explained later on, at room tem- 
perature radiative processes dominate over nonradiative 
processes and therefore, for the study of nonradiative pro- 
cesses, it is necessary to measure the PL decay at low tem- 
peratures. Typical PL decay curves, measured for H-PSi at 
a photon energy of 2.03 eV (610 nm) and at various tem- 
peratures, are presented in Figure 2. A pronounced de- 
pendence of the PL decay on temperature can clearly be 
seen, similar to the results of other groups [1,2,53]. As the 
temperature decreases, the PL decay time becomes signifi- 
cantly longer (by two orders of magnitude over the entire 
range of measured temperatures). Notice that the tem- 
poral behavior of the PL cannot be described by a simple 
exponential decay function (see the semi-logarithmic scale 
of Figure 2) and is typically fitted to a stretched exponen- 
tial decay function [54,55]. This nonexponential decay is 
common to disordered systems and has been attributed to 
a dispersive diffusion of the photo-excited carriers [54]. 
The solid lines in Figure 2 represent the best fit of the PL 
decay curves to a stretched exponential function, given by 

I(t)=I 0 exp[-(t/rf] (1) 

where r is the PL lifetime, and ft is the dispersion expo- 
nent that was found to vary in between 0.4 to 0.8 and 
will not be discussed here (see [37] for more details). 
Arrhenius plot (semi-logarithmic scale versus the in- 
verse temperature) of the measured PL lifetime for both 
H- and O-PSi (at a photon energy of 2.03 eV) is shown 
in Figure 3a, presenting exponentially fast decays at 
high temperatures and approximately long and constant 
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Figure 3 PL lifetime and integrated PL (a) Arrhenius plot of the 
PL lifetime (on a semi-logarithmic scale) as a function of 1/7", at a 
photon energy of 2.03 eV (610 nm) for H-PSi (red circles) and O-PSi 
(black squares). The solid lines represent the best fit to the singlet- 
triplet model of Eq.2. (b) Arrhenius plot of the integrated PL. Inset 
shows the schematics of the excitonic two-level model with the 
upper excitonic singlet-triplet state and the ground (no exciton) 
state. The arrows represent the allowed (from the singlet) and the 

forbidden (from the triplet) optical transitions, 
k J 

decay times at low temperatures. This unique behavior 
of the PL decay has been attributed to a splitting of the 
excitonic ground state (i.e., the photo-excited electron- 
hole pair) due to the Coulomb exchange interaction, 
giving rise to a lower triplet level (S = 1) and an upper 
singlet level (S = 0) [53] (see inset to Figure 3b). The 
splitting is further enhanced by confinement of the 
charged carries in small nanocrystals, giving rise to a 
larger excitonic overlap. Optical transitions from the 
lower triplet and the upper singlet states are forbidden 
and allowed respectively, due to spin selection rules 
[1,2,39,40,53]. However, the lifetime of the triplet state 
becomes weakly allowed due to spin-orbit interaction 
[39,40,53]. Hence, the triplet lifetime is expected to be 
considerably longer than the singlet lifetime. At low 
temperatures (where kT < < A, and A is the singlet- 
triplet splitting energy; see inset to Figure 3b), only the 
triplet level is populated and therefore, the PL decay 
time is dominated by the triplet lifetime and is relatively 
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long (the low- temperature plateau regions in Figure 3a). 
As the temperature increases (above approximately 
30 K), the upper singlet level becomes thermally popu- 
lated and the overall lifetime shortens according to the 
following expression: 



1 = f i+ ±exp(-A/kT) 
t r g+ exp(-A/£r) 



(2) 



where r R stands for the radiative decay time and r L and tjj 
are the lower triplet and the upper singlet lifetimes respect- 
ively (g = 3 is the levels degeneracy ratio) [1,39,40,53]. At 
high temperatures, the decay time is dominated by the 
much faster upper singlet lifetime. 

From Figure 3a we found that within the experimental 
errors, the upper singlet decay times of H- and O- PSi 
(at photon energy of 2.03 eV) are essentially the same 
(1.0 ± 0.2 us and 1.3 ± 0.2 us for H-PSi and O-PSi, re- 
spectively). However, at low temperatures the H-PSi 
decay time is faster than that of the O-PSi (200 ± 50 us 
relative to 480 ± 50 us, respectively). To further explore 
the differences between H- and O- PSi decay times, the 
singlet and the triplet lifetimes as well as the energy 
splitting were extracted over the measurement's range 
of photon energies and are plotted in Figure 4. As ex- 
pected, the upper singlet lifetime (ru) is significantly 
shorter (by about one to two orders of magnitude) than 
the lower triplet lifetime (r L ) over this range of photon 
energies. The energy splitting, AE (see Figure 4b), was 
found to vary between 7 and 25 meV, as in [53] and in 
accordance with the calculated values in [56-58]. Com- 
paring H- and O- PSi, we note that the upper singlet 
lifetimes and the excitonic energy splitting of both H- 
PSi and O-PSi remarkably coincide over the entire range 
of measured photon energies (see Figure 4a,b), while 
the lower triplet lifetime of H-PSi is shorter than that 
of O-PSi over the same range of energies (Figure 4c). 
This result is the basis for our conclusion (to be dis- 
cussed hereafter) that oxidation of (freshly prepared) 
H-PSi gives rise to slower nonradiative lifetimes, leaving 
radiative lifetimes unaffected. 

Discussion 

As explained above, the main finding of this work is that 
the oxidation of freshly prepared luminescent PSi gives 
rise to slower triplet lifetimes, keeping the upper singlet 
lifetimes unaffected. Before discussing the implications 
of this result, let us denote that the measured decay rate 
is the sum of two competing relaxation processes given 
by 
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Figure 4 Triplet and singlet lifetimes and energy splitting, (a) 

the upper singlet lifetime; (b) the excitonic energy splitting; (c) the 
lower triplet lifetime (extracted from the fit to the singlet-triplet 
model; see Figure 3) as a function of the photon energy. 



where Tr 1 is the radiative transition rate (given by Equation 2), 
Tnr is the nonradiative relaxation rate, and r _1 is the total 



decay rate. The integrated PL (i.e., the area below the PL 
spectrum shown at the inset to Figure 1) is proportional to 
the quantum yield that is given by the ratio of the radiative 
to the total decay rate, rj = r^/r -1 = t/tr. The variation 
of the integrated PL with temperature is shown in Figure 3b 
on a semi-logarithmic scale, similar to that of Figure 3a for 
the PL lifetime. Notice that while the PL lifetime varies by 
approximately two orders of magnitude over the 30 to 
300 K temperature range, the integrated PL varies by less 
than 3. Hence, one concludes that at this temperature 
range, r R < < r NR , leading to, r ~ r R (Equation 3), and rj ~ 
constant (as in reference [37]). Thus, at temperatures above 
30 to 40 K the measured lifetime is dominated by radiative 
transitions. In addition, the strong dependence of the upper 
singlet lifetime on photon energy (a decrease from 6 to 7 (is 
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at 1.6 eV down to 200 to 300 ns at 2.3 eV; see Figure 4a), 
suggests again that this lifetime should be associated with 
radiative transitions (where r u ~ r R u < < r NR u ). In this 
case, the fast radiative lifetime is due to the influence of 
confinement on the spontaneous emission rates in small 
Si nanocrystals [39,40]. On the other hand, the lower trip- 
let lifetime that is dominant at low temperatures is ap- 
proximately constant (varies by less than factor of 2 over 
the same range of energies) and roughly independent of 
the photon energy that probes a given size of nanocrystals. 
This suggests that the low-temperature lifetime should be 
associated with a nonradiative relaxation time (of the 
whole system) that dominates over the (forbidden) ra- 
diative triplet lifetime, in agreement with [37,39,40]. 

Turning back to our main findings, we conclude that 
oxidation (in ambient conditions) has a minor impact 
on the size of the nanocrystals (giving rise to about 3% 
blue shift of the PL spectrum) and no noticeable effect 
on the radiative lifetime and the excitonic energy split- 
ting (via their dependence on photon energy). On the 
other hand, nonradiative relaxation times, which are as- 
sociated with the state of the surface, are expected to be 
sensitive to oxidation and to a modification of surface 
bonds as experimentally observed (see Figure 4c). This 
result can be explained by the EV model [39,40], which 
assigns the slow nonradiative relaxation times to reson- 
ant coupling between surface vibrations and quantized 
electronic sublevels in the conductance/valence bands 
of the nanocrystals. The stronger is the coupling be- 
tween these electronic states and surface vibrations, the 
slower are the nonradiative lifetimes [39,40]. Hence, ac- 
cording to this model, the longer lifetime measured for 
O-PSi (compared to H-PSi) should be assigned with the 
larger electronegativity of oxygen (relative to hydrogen) 
that gives rise to larger dipole strength of the Si-O-Si vi- 
bration [47]. 

Finally, let us point out that the conclusion r R < < r NR 
(for both types of PSi; see Figure 4) implies that the PL 
quantum yield should approximately be constant. This 
conclusion provides a simple explanation to the slight 
variation of the PL intensity under oxidation, as oxida- 
tion modifies nonradiative relaxation times associated 
with the PSi surface. However, this has a minor impact 
on the PL quantum yield as the PL efficiency is practic- 
ally independent of the nonradiative relaxation times at 
high temperatures [39,56,57] and is mostly influenced by 
the nanocrystals size [59,60]. 

Conclusions 

In conclusion, using temperature-dependent, time-resolved 
PL spectroscopy for probing both radiative and nonradia- 
tive relaxation processes in freshly prepared and oxidized 
PSi, we were able to show that radiative processes should 
be associated with quantum confinement in the core of 



the Si nanocrystallites and therefore, are not affected by 
oxidation. On the other hand, nonradiative relaxation pro- 
cesses are affected by oxidation and by the state of the 
nanocrystallites surface. These results are consistent with 
the extended vibron model that assigns radiative relax- 
ation to QC, while nonradiative processes are assigned to 
surface chemistry. 
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